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Abstract 



Neutrino induced reactions on C, an ingredient of liquid scintillators, 
have been studied in several experiments. We show that for currently avail- 
able neutrino energies, E v < 300 MeV, calculated exclusive cross sections 
12 C ss (z/, l) 12 N gs for both muon and electron neutrinos are essentially model 
independent, provided the calculations simultaneously describe the rates of 
several other reactions involving the same states or their isobar analogs. The 
calculations agree well with the measured cross sections, which can be there- 
fore used to check the normalization of the incident neutrino spectrum and 
the efficiency of the detector. 

PACS numbers: 13.15.Dk, 21.60.-n,23.40.Bw,25.30.Pt 
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A solid theoretical understanding of cross sections in neutrino-induced reactions on light 
nuclei is becoming a necessity, in particular for 12 C, an ingredient of liquid scintillators, 
and for 16 0, the basic component of water Cerenkov detectors. Detectors containing these 
isotopes are used to measure fluxes of atmospheric and supernova neutrinos, and in searches 
for neutrino oscillations. These tasks to a large extent require knowledge of the corresponding 
cross sections, which often have not been measured. When experimental cross sections are 
available it is therefore important to compare them with calculations. 

For 12 C a number of experimental results exist. These include measurements of charged- 
current reactions induced by both electron- [T|] and muon- neutrinos [0-f|], exciting both the 
ground and continuum states in 12 N. The inclusive cross section for 12 C(y e , e) 12 N*, measured 
by Karmen Q, LAMPF (with large errors) ||, and recently by the LSND collaboration [|J, 
agree well with calculations. [A previously reported measurement of both exclusive and 
inclusive neutrino cross sections on 12 C [@ appears to be inconsistently large 0; we will 
not consider that measurement here.] By contrast, there is a disturbing discrepancy be- 
tween calculations and the LSND value of the cross section for the inclusive reaction 
12 C(z/ M , /i) 12 N*, which uses higher energy neutrinos from pion-decay in flight []3]|||| . The dis- 
agreement is disturbing^ in light of the apparent simplicity of the reaction and in view of the 
fact that parameter-free calculations, such as those in |)],|7j describe well other weak processes 
which are governed by the same weak-current nuclear matrix elements. Here we examine 
whether similar problems affect our understanding of the exclusive reactions 12 C(u e , e) 12 N 9 <, 
IHII and 12 C(z/ A1 , /x) 12 N 5S 0,f|||. The exclusive process is a useful monitor of the neutrino 
flux. If calculations are reliable and reproduce experiment, then the normalization of the 
experimental neutrino flux must have been correctly modeled. 

A calculation of the exclusive cross section can be tested by computing rates of related 
processes and comparing to data. The 1 + T = 1 ground state in 12 N is the analog of the 15.11 
MeV state in 12 C and of the ground state of 12 B. This allows us to use (3 + decay from 12 N 
back to the 0+ T = ground state of 12 C, muon capture from 12 C to the ground state of 12 B, 
ineleastic electron scattering to the 15.11 MeV state in 12 C, and Ml decay from that state 
to the ground state to calibrate elements of the calculation. In what follows we calculate 
the exclusive neutrino cross sections in several ways to see whether when all the above data 
are reproduced, the different models can produce significantly different predictions. 

Our first approach is a restricted shell-model calculation. Describing all the above reac- 
tions is straightforward in the shell model as long as we can neglect contributions of two-body 
operators (i.e. meson exchange). The reduced matrix element of an arbitrary operator O is 
then given by 

(i + A\\d\\o\o) = Y,(j\\d\\f)OBD(j,f) , (i) 

hi' 

where the one-body transition densities are defined by 

OBD(j,f) = (l + ,l||[a^-] J ^ T = 1+ ' 1 ||0 + ,0) . (2) 



1 The disagreement is not universal. In [1C] a cross section that agrees with LSND is obtained 
using the elementary particle method. We discuss the applicability of the approximations used in 
[10] elsewhere [11]. 



2 



If we further assume that all structure in the low-lying states with J 71 = 1 + ,T = 1 is 
generated by the eight "valence" nucleons moving in the p-shell, there are only 4 one-body 
densities and 4 single-particle matrix elements (j| |0| Ij'), which contain all the momentum- 
transfer dependence and that are simple to calculate. Futhermore, the Gamow- Teller (GT) 
matrix element for the f3 + decay of 12 N, the Ml radiative width of the state at 15.11 MeV 
in 12 C, and the form factor for the excitation of this state in electron scattering jl2| depend 
only on three one-body densities; they are independent of the combination OBD(p 3 / 2 ,pi/2) + 
OBD(pi/2iPz/2). (See the review |R| for a general discussion.) 

The most straightforward way of obtaining the one-body densities is by diagonalizing a 
thoroughly tested residual interaction, such as the one given in Ref. [14"f. It is well known, 
however, that the resulting p-shell one-body densities do not describe the reactions above 
very well. Typically, the calculated GT strength is too large (the origin of "GT quenching") 
and the electron scattering form factor is too high in the first lobe. In a number of papers 
P~5|— rS| over the last twenty plus years, however, it has been shown that one can modify 
the one-body densities (ad hoc) in such a way that the three experiments are correctly 
reproduced, with the form factor adequately described up to the first minimum at |g| ~ 1.5 
fm -1 . This is the approach we follow here, adjusting the densities to reproduce the data. The 
so-far undetermined fourth combination of densities can be fixed (with some uncertainty) 
by the muon capture rate to the ground state of 12 B [13]. This process is the only one that 
tests the momentum dependence of the axial current, since it takes place at q 2 = —0.74m 2 . 
For our analysis we take an averaged experimental rate, u>(l + ) = 6200 ± 200 s _1 [|19j, and 
use the Goldberger-Treiman relation for the induced-pseudoscalar coupling constant. 

Line 2 of Table I contains the resulting one-body densities, adjusted here to reproduce all 
the data discussed above. Lines 3-5 contain one-body densities used by other authors, and 
constrained to different subsets of the above data. To achieve an overall agreement with all 
the data, we renormalize, in addition, either the weak axial form factor Fa or the magnetic 
form factor Fm- Haxton |T5[] (line 3) required that the j3~ decay of 12 B, rather than 12 N 
(but not the muon capture), be well described; because of isospin violation in the ft values, 
his densities require a renormalization factor of 0.873 for the axial current form factor Fa to 
fit the 12 N decay. The entries in line 4 are based on the extreme single particle model |16 



in which 12 C is represented as a closed p 3 / 2 subshell, and 12 N or 12 B have one nucleon in pi/ 2 
subshell. Here a renormalization of 0.414 is required for Fa, and 0.484 for the magnetic form 
factor Fm- Finally, line 5 contains the one-body densities of Ref. [H]], based on a minimally 
modified Cohen-Kurath interaction. They require only a small Fa renormalization of 0.925. 
In all cases harmonic oscillator wave functions are used, with slightly different values of the 
oscillator parameter b taken from the original references. The table clearly illustrates that 
the extraction of one-body densities from data is not a unique procedure; it depends on 
other assumptions about wave functions, etc. We discuss the effects of these differences on 
the exclusive neutrino cross sections shortly. 

While two-body effects, i.e. meson exchange corrections, are expected to be relatively 
small (5-10 %) for the momentum transfer considered here [GO], configurations beyond the p 
shell might explain the need for a drastic renormalization of the one-body densities produced 
by a reasonable p-shell Hamiltonian. We therefore calculate the rates of all the reactions 
above, including exclusive neutrino capture, in the Random Phase Approximation (RPA), 
which does include multishell correlations, while treating the configuration mixing within 
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the p shell only crudely. The same method, extended to the continuum, has been used 
in the calculation of the neutrino charged and neutral current inclusive scattering [|6j|7l,|2ll , 
(disagreeing disturbingly with one experiment, as noted above) and muon capture [fZ2] . 



To compensate for the crude description of the p-shell dynamics we use an overall quench- 
ing factor of 0.258 by which we multiply the rates, respectively the cross sections of all 
processes under consideration, for all momentum transfers. With this multiplicative factor 
the fl + decay, Ml width, electron scattering form factor and partial muon capture rate are 
all adequately described, and the exclusive cross sections, discussed below, are readily cal- 
culated. (It has been known for some time that the RPA is capable of describing the shape 
of the (e, e') form factor for the 15.11 MeV state ) 

Our third approach is the "elementary-particle treatment" (EPT). Instead of describing 
nuclei in terms of nucleons or other constituents, the EPT considers them elementary and 
describes transition matrix elements in terms of nuclear form factors deduced from experi- 
mental data and constrained by transformation properties. The EPT was implemented in 
the A = 12 system in Refs. [p4Tp6| for neutrino energies up to E v = 100 MeV. Here we 
extend the approach to the higher neutrino energies relevant to the LSND decay-in-flight 
z/^'s by appropriately including the lepton mass, which was ignored in Ref. ||24j| , in the kine- 



matics. A nonzero lepton mass requires, in turn, that the induced pseudoscalar term, also 
neglected in Ref. |24j], must be included as well. Here we used a modified form of Eq. (16) 
in Ref. kindly furnished to us by Professor Kubodera and for the F P {Q 2 ) we used the 
simple ansatz: 
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Fp(Q 2 ) = -^*U(0) 



where, as usual, Q 2 = —q 2 , and the empirical parameter £ has been fixed from muon capture. 

We turn now to the evaluation of the exclusive cross sections in the three approaches. 
Within the shell model and the RPA the cross section is easily evaluated once the one-body 
densities and free-nucleon form factors (for which we use standard values) are specified; the 
general prescription can be found, e.g., in Ref. |27|| . In the EPT the evaluation is even 
simpler since the nucleus is elementary. One complication often not considered is related to 
the Coulomb interaction of the outgoing charged lepton. The usual treatment, as in nuclear 
(3 decay, involves the Fermi function F(Z, E), which is the ratio of the Coulomb continuum 
s-wave and the corresponding free s-wave wave. This approximation is valid, however, only 
for lepton momenta pR < 1, where R is the nuclear radius. For 12 C it can be used for 
electrons with up to about 50 MeV of kinetic energy, but it is justified only for muons of 10 
MeV or less of kinetic energy. As energy increases, F(Z,E) approaches a constant value, ~ 
1.17 (Z = 7) for v reactions on 12 C, and ~ 0.90 (Z = 5) for //-induced reactions on 12 C. On 
the other hand, it is intuitively clear that as the lepton energy becomes much larger than 
the Coulomb energy, the Coulomb correction should approach unity. In order to keep things 
simple, and since the correction is a relatively minor one, we scale at higher energies the 
lepton cross section by the ratio 

p eff (E+(V)) , „ (V), /t r\ .Mat 

pE > where Peff = Pi 1 + V } ' {V) = ± ^2R ' ( } 

and (V) represents the average Coulomb potential. The two approximations are smoothly 
connected. 
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The results appear in Table II. For the ^-induced reaction with neutrinos from the 
muon decay-at-rest |lj the agreement between the experimental and calculated exclusive 
cross section is perfect in all the models. This cross section, corresponding to an average 
momentum transfer of only about 50 MeV, appears totally fixed by the requirements we 
impose on each of the calculations. For the ^-induced reaction the average momentum 
transfer is about 150 MeV. But even in this case the different variants of the shell model, 
the RPA and the EPT give quite similar cross sections. Again, the agreement with the 
LSND experimental value is good. 

In Figures 1 and 2 we show the exclusive cross sections for 12 C(V e , e) 12 N 9S and 
12 C(z/ M , /i) 12 N 9S , respectively, as a function of neutrino energy. For z/ e 's the cross section 
is essentially the same in all the approaches we consider. The agreement between the vari- 
ous models is also quite good for z/^'s. Moreover, after a very steep rise from the threshold, 
this cross section quickly saturates and becomes more or less independent on the neutrino 
energy. Thus, the exclusive cross section is simply proportional to the total number of u^s 
above threshold. However, it is essentially independent of their energy distribution. 

For completeness, Table III lists the cross sections for neutral-current excitation of the 
15.11 MeV, 1 + ,T = 1 state in 12 C, which was measured by Karmen jl| for neutrinos from 
decay at rest. Once again all the models agree with one another and the data. The neutral- 
current cross section for v^'s from decay-in-flight and for i? M 's are probably not measurable, 
but unsurprisingly the calculations continue to agree with one another. 

In conclusion, we have shown that the exclusive cross sections for neutrino energies 
available now are calculable in a variety of ways, with results that are nearly identical, 
and agree very well with the data. These processes can therefore be used as a check of 
the neutrino flux normalization, and detector efficiency and indicate that the discrepancy 
between the measured and calculated inclusive 12 C(z/^, /z) 12 N* cross sections is not caused 
by an underestimate of the neutrino flux normalization. The source the disagreement must 
be found elsewhere. 

We would like to thank Professor Kubodera for his help with the EPT. We also appreciate 
numerous discussions with the members of LSND collaboration. We were supported in part 
by the U.S. Department of Energy under grants DE-FG05-94ER40827 and DE-FG03-88ER- 
40397, by the U.S. National Science Foundation under grants PHY94-12818 and PHY94- 
20470, and by the Swiss Nationalfonds. 
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TABLES 



TABLE I. Fitted one-body densities OBD(j,f). The parameter b is the oscillator length. The 
range in line 2, column 4, reflects the uncertainty in the muon capture rate; in other lines only the 
most probable capture rate was used. 



J, J 



1/2,1/2 



1/2,3/2-3/2,1/2 



1/2,3/2+3/2,1/2 



3/2,3/2 



b (fm) 



present 
Ref. |5| 
Ref. |6| 
Ref. |7| 



-0.113 
-0.111 
0.0 
-0.27 



0.106 
0.337 
1.0 
0.318 



0.666±0.4 
0.875 
1.0 
1.03 



0.24 
0.086 
0.0 
0.116 



1.67 
1.76 
1.77 
1.82 



TABLE II. Comparison of calculated and measured cross sections, in units of 10 cm 2 
and averaged over the corresponding neutrino spectra, for the neutrino induced transitions 
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c 



► 12 N flfl and 12 C gs 



,12 



B ffS . For the decay at rest the v e spectrum is normalized from 
E v = while for the decay in flight the Vn and spectra are normalized from the corresponding 
threshold. The cross section for v e is not quoted since the decay-at-rest neutrino source does not 
contain any P e 's. 



12 C(z, e ,e-) 12 N, s 


12 C( W -) 12 N SS 


12 C(^,/i+) 12 B 5S 


decay at rest 


decay in flight 


decay in flight 



experiment 
experiment 
OBD of this work 
OBD of Ref. Ill 



OBD of Ref. gH 
OBD of Ref 
CRPA @| 
EPTtt JH 



8.9±0.6±0.75 
9.1±0.4±0.9 [ 
9.1 ± 0.1 

8.8 

9.4 

9.4 

8.9 

9.2 



64±io ± io m 




63.5 1 " ±5 


24 


60.4 


23 


65.4 


22.6 


62.3 


23.7 


63.0 


26 


62.9 


21.5 



t The uncertainty reflects the range corresponding to the uncertainty in the muon capture 
rate. A similar range is presumably valid for the other approaches in lines 4-6. 
ft Extended to muon neutrinos, K. Kubodera, private communication. 
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TABLE III. Comparison of calculated and measured cross sections, in units of 10 _42 cm~ 2 
and averaged over the corresponding neutrino spectra, for the neutrino-induced neutral-current 
transitions 12 C gs (u, z/) 12 C (15.11). In column 2 the contributions from v e and Un are added. All 
spectra are normalized "per neutrino" i.e. from E u =0. 



12 C(^/^eK) 12 C(15.11) 12 C(^,^) 12 C (15.11) 12 C(^) 12 C (15.11) 



decay at rest 



decay in flight 



decay in flight 



experiment 
OBD of this work 



OBD of Ref. (TJ 
OBD of Ref. g 
OBD of Ref |T 
CRPA M 



EPT H 



11±1.0±0.9 
9.3 
9.3 
9.4 
9.6 
10.5 
9.4 



22.3 


15.2 


21.9 


15.3 


23.7 


16.4 


22.9 


15.9 


27.2 


17.5 


24.3 


14.5 
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FIGURES 



FIG. 1. Cross section for 12 C(f e , e) 12 N gs as a function of the v e energy, in units of 10 2 cm 2 . 
The full line is for the shell model with one-body densities from line 2, Table I, the dashed line is 
the modified Cohen-Kurath one body densities from line 5, Table I JlTj (the curves for the other 
shell model variants are very similar and are not shown), the dot-and-dashed line is for the RPA 
0, and the dotted line is for the EPT. 

FIG. 2. Cross section for 12 Cz^, n) 12 N gs as a function of the i/„ energy, in units of 10 _42 cm 2 . 
The notation is as in Fig. p]. 
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